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Abstract. We present a finite-dimensional version of the quantum model for the stock 
market proposed in [C. Zhang and L. Huang, A quantum model for the stock market, 
Physica A 389 (2010) 5769]. Our approach is an attempt to make this model consistent 
with the discrete nature of the stock price and is based on the mathematical formalism 
used in the case of the quantum systems with finite-dimensional Hilbert space. The rate 
of return is a discrete variable corresponding to the coordinate in the case of quantum 
systems, and the operator of the conjugate variable describing the trend of the stock 
return is defined in terms of the finite Fourier transform. The stock return in equilibrium 
is described by a finite Gaussian function, and the time evolution of the stock price, 
directly related to the rate of return, is obtained by numerically solving a Schrodinger 
type equation. 



1 Introduction 

The methods of quantum mechanics are among the tools used in order to get a deeper 
insight in the complexity of the financial markets [1-9] . The time evolution of the price 
and the trend of the stock in the financial markets depends on many factors including 
the political environment, market information, economic policies of the government, 
psychology of traders, etc. Prices at which traders are willing to buy or sell a financial 
asset are not more determined by the development of industry, trade, services, situation 
at the market of natural resources and so on [7] . The information exchange and market 
psychology play important role in price dynamics. 

The mathematical modeling of price dynamics of the financial market is a very com- 
plex problem. We could never take into account all economic and non-economic con- 
ditions that have influences to the market. Therefore, we usually consider some very 
simplified and idealized models, a kind of toy models which mimic certain features of a 
real stock market. For a better understanding of the complexity of the financial markets 
one has to use several complementary models. 

The variables describing the market dynamics are discrete variables, not continuous 
ones. For example, the stock price at a given moment of time is an integer multiple of 
a certain minimal quantity, a sort of quantum of cash (usually, 1/100 or 1/1000 of the 
currency unit) [5]. The time evolution of the price is thus described by a step function. 
The quantum models facilitate the description of phenomena not fully explained by the 
classical theory of prices. The lack of simultaneously observability of certain variables 
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and the interference of attempted values noted by traders and economists represent such 
quantum effects. There is the impossibility of observing prices and their instantaneous 
forward time derivative even if we consider price as a continuous variable [10) . The rate 
of return being directly related to the price, there is also the impossibility of observing 
the rate of return and its instantaneous forward time derivative. The trade of a stock 
can be regarded as the basic process that measures its momentary price. The stock price 
can only be determined at the time of sale when it is between traders. We can never 
simultaneously know both the ownership of a stock and its price [HI [H] • The stock price 
p is a variable taking numerical values but the ownership seems to be a very different 
kind of variable. 

Our aim is to present a quantum model describing the time evolution of the rate of 
return. We shall take into consideration only the case of stock markets with a price limit 
rule: the rate of return 

Po Po 

in a trading day can not be more than ±q% comparing with the previous day's closing 
price po, that is, we have 

100 ~ ~ 100 

For example [5], in most Chinese stock markets q — 10. In their activity, the traders 
do not take into consideration an arbitrary small rate of return. They approximate the 
rate of return by integer multiples of a minimal one. We shall consider 1Z as a discrete 
variable with the only possible values 

2_ _?zl 3_ n -J— £i _2_ 

100' 100' '" ' 100' ' 100' '" ' 100' 100 

and by following the analogy with quantum mechanics, we describe the rate of return at 
a fixed moment of time by a wave function 

■ I ioo ' ioo ' ••■ ' 100 ' 100 J ' ^ 
such that |$ ( j^q) 1 2 is the probability to have a return rate equal to ^g. 

2 Infinite-dimensional quantum system 

In the case of a quantum particle moving along a straight line, the possible positions form 
the set R=(— oo, oo), and each state of the system is described by a function ifj-.M. — >C 
satisfying the relation (and called square integrable) 

|^(x)| 2 dx < oo. 

> 

The space of square integrable functions considered with the scalar product 



(ipWh) = ipi (x) dx (2) 

J — OO 

is an infinite dimensional Hilbcrt space, denoted by L 2 (M). To each function — 
we associate its Fourier transform Thp] :R — denned as 



F[tp](p) = -= e~ ipx/h ij(x)dx (3) 
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where H is Planck's constant h divided by 2w. The transformation tp t— > T[tp] is a unitary 
transformation 

(Jfyi], F[ip 2 }) = (lh,ih) for any ipi,ip2 € L 2 (R) 
and its inverse is the transformation 

1 f°° 

F-\4]{x) = -== e^l h ^{p)dp. (4) 
V^7rft j-oo 

The normalized function corresponding to a function tp ^ 0, namely, 



satisfies the relation 



*:M — >C, #(x) = — tb(x) 



/oo 
\*(x)\ 2 dx = 1. 
-co 



In the state described by the normalized function 'J, the number 

rb 

\^(x)\ 2 dx 

represents the probability to find the particle in the interval [a, b], and 

\nmp)\ 2 dp 



represents the probability to have a momentum in the interval [a,b]. 

Each observable A is described by a linear operator A satisfying the relation 

(Ai> 1 ,,k) = (ii> 1 ,Aih) (5) 

called a Hermitian operator. The position is described by the operator 

ipt-txip where (xip)(x) = xip{%) (6) 
and the momentum by the operator 

ip i — y p tp where pip = —\h-^-. (7) 

ax 

The position and momentum operators satisfy the relation 

p = T~ x xT. (8) 

In a state described by the normalized wavefunction "J, the numbers 

(x) = {if!, x and (p) = 

represent the mean value of coordinate and momentum, respectively. The time evolution 
of the state of the system satisfies the Schrodiger equation 

i&*{x,t)=(±p + V(x,tj)*( Xi t) (9) 

where \x is the mass of the particle and the operator V(x, t) is obtained from the potential 
energy V(x,t) by using the substitution x i— > x. 



3 Finite-dimensional quantum systems 



We obtain a very simplified version of the infinite-dimensional quantum system presented 
in the previous section by assuming that we can distinguish only a finite number d of 
positions for our particle. It is convenient to consider that rf = 2s + l is an odd number 
and to assume that 

_ ~ ~ ~ (10) 



is the set of all the positions we can distinguish. From the relations 



lim 

d— yoo 



= 0. 



lim {—s) 

d— voo 



d— >oo 



it follows that, in a certain sense, we have lim Sd — (-co, oo). 

d— > oo 

Each state of our system is described by a function i\)\S d — 
these functions considered with the scalar product 

(^1,^2) = E ^l( x ) fo( x ) 



The space T~L of all 
(11) 



is a d-dimensional Hilbert space isomorphic to C d . To each function ijj : Sd 
associate its Fourier transform F[ijj] :Sd — >C defined as [12 ) [T3 ] 

x£S d 

that is, 

m (k^f) = ^ t_f^ kn 1> ■ 
The transformation tp 1— > F[tp] is a unitary transformation 

(F^.FM) = {ipufa) for any ipufoeH 
and its inverse is the transformation 

The normalized function corresponding to a function -0 7^ 0, namely, 



we 



(12) 



*(x) = 



satisfies the relation 



1 — — s * * 



= 1. 



(13) 

(14) 
(15) 



In the state described by the normalized function 'J, the number 



represents the probability to find the particle at point bJj, and the number 



F[*]lkJ% 
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the probability to have a momentum equal to fcy ^f-- Each observable A is described by 
a linear operator A satisfying the relation 

(Ml, ^2) = (ipi,Aip 2 ) for any i/) b i/j 2 eH 

called a Hermitian operator. The position is described by the operator 

Q :H — > U : ijj ^ Qijj 

where 



(Qip)(x) = xtp(x) 



that is, 



(16) 



The function -0 is a function of discrete variable and the expression ^ is meaningless. 
Therefore, we define the operator corresponding to the momentum by following the 
analogy with the relation as 



P-.n^n, P = f- x qf. 

In a state described by the normalized wavefunction "J, the numbers 



(17) 



(Q) = (%QH()= £ nJ% 



and 



(P) = (*,P*) = (^P^QP*) = (P[*],QP[*]) 



represent the mean value of coordinate and momentum, respectively. The time evolution 
of the state of the system satisfies a Schrodiger type equation 



where the operator 



d 

i — f = P*. 

dt 



H = ±P> + V(Q,t) 



(18) 



(19) 



is the corresponding Hamiltonian. If H does not depend on time than the solution of 
the Schrodingcr equation satisfies the relation 



* [n«/^,t] =e 



(20) 



The finite-dimensional quantum system presented in this section converges for d — > 00 
to the system presented in the previous section, in a certain sense. 
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4 A quantum model for the stock market 

We consider a stock market with a price limit rule [8] , for which the rate of return 1Z in 
a trading day can not be more than ±q%. By assuming that 

r, = {^± 2nk _2_\ foil 

I 100' 100 ' ■•• ' 100' 100 / K^^-J 

is the set of all the possible values of the rate of return we use a function ip : Cd — > C 
in order to describe the rate of return. The space K, of all the functions (p : Cd — > C 
considered with the scalar product 

8 

I 

n=—q 



(<pM = E Mito)Mito) (22) 



is a Hilbert space of dimension d = 2q+l, isomorphic to C d . To each function tp-.Cd^C 
we associate its Fourier transform F[ip] :Cd — >C defined as 

^Hi§o) = 72 £ " fo XiM)- (23) 

n= — q 

The transformation <p H ► .F[y] is a unitary transformation 

(F[^i],F[^ 2 ]) = (<Pi,<P2) for any <£>i,<p 2 €£ 
and its inverse is the transformation 

F - 1 M(m) = 73 £ e¥t >(M)- (24) 

fc— — q 



The normalized function corresponding to a function ^ 0, namely, 

1 



$:£ d ^C, $(^)= 1 ^(-n-) 



satisfies the relation 

E I^(t§o)| 2 = 1 - (25) 

n— — q 

By following the analogy with quantum mechanics, we assume that in the case of a rate 
of return described by the normalized function <£>, the number 

I ^ v 100 y I 

represents the probability to have a return rate equal to . 

In our model the financial variables are described by Hermitian operators. The rate of 
return corresponds to the position, and by following the analogy with quantum mechanics 
we define the corresponding operator as 

Mean value of the rate of return 1Z in the case of a normalized function $ is 

<i?) = <M$)= E I*(im)| 2 - 

n— — q 

The rate of return has been defined by using the stock price as 

n = P~ Po = P_ _ 1 

Po Po 
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This relation, written in the form 

p = po + po n 

allows us to define the price operator as 

p = p l + po R (27) 

where I : K, — > IC, lip = ip is the identity operator. 

In our discrete version, we cannot define the operator T corresponding to the trend of 
return rate by using the method from [8 as a differential operator. We define it directly 
by using the finite Fourier transform as 

f-.K. — >K, f = F^RF. (28) 

The set of functions {5 n }1 l= _ ql where 

is an orthonormal basis in JC, called the canonical basis. The function 5 n is an eigen- 
function of the operator R corresponding to the eigenvalue 

RS n = m 5 ^ ( 3 °) 

and any function ip : Cd — > C can be written as a linear combination of S n 

<P= i Hm) 5 n- (31) 



The set of functions {^ n }n=- q defined as = F [S n ], that is, 



as 



$„ : C d — > C, $„ (4) = ^e^ fc " (32) 

is a second orthonormal basis in K,. Each $„ is an eigenfunction of T 

f$ n =F- 1 RFF- 1 [6 n ]=F- 1 RS n = F _1 [<5„] = ^ $„ (33) 

and any function ip : Cd — > C can be written as a linear combination of $„ 

<? 

n— — q 

Since 

($ B)V ) = {F-\8 n U) = (8 n ,F[<p\) = F[ip] (^) 
each function ip : Cd — > C admits the representation 

v= £ F l<p](m) *»• (34) 

n— — g 

In the case of a normalized function $ we have 

i \m(M 2 = t i(^,*>i 2 = i 

k——q k——q 
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■7o. 2(2:) 



7l (a:) 



72(a) 



Figure 1: The functions 70.2 (left), 71 (center) and 72 (right) in the case q = 10. 



and the number 



!*"[*] (155) I 



is the probability to have a trend equal to ' 



100 • 



By following the analogy with the quantum systems, we assume that the function $ 
describing the rate of return satisfies a Schrodinger type equation 



where the operator 



H = 7-P 2 + V(Q,t) 
2/i 



(35) 
(36) 



is the corresponding Hamiltonian. By denoting 

4> n = {5 n ,§) and h n>m = (S n ,HS m ) 
the Schrodinger equation can be written in the matrix form 



dt 



( 4>- s {t) \ 



h-s+l-s fr-s+l,-s+l 



h- s ,s 
h-s+i 



+l,s 



V h. 



( 4>-s{t) \ 

<t>-s+l{t) 



that is, it is a system of d first order differential equations with d unknown functions 
<f>-a, (f>-s+i, ■•■ , (f>s- 



5 Market information and the price evolution 

The rate of return of the stock market in equilibrium is usually described by a Gaussian 
function [8] [14] . In the case of our model we use a function 

7 Q : L d — >R, 7a(x) = 1 g a {x) (37) 

V\9a,9a) 

(see Figured]) depending on a parameter a£ (0, 00), where 

9a (tm) = £™=-oc e^(™ d +«) 2 = ^ 3 (5, • (38) 
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By using the properties of the Jacobi function 



6 3 (z,t) = e™ Ta ~ e 27riaz (39) 

a.—— oo 

Ruzzi [15] has obtained the relation 

F\g a ] = ^=gi. (40) 

In the case of a Hamiltonian 

H=±f 2 + V(R,t) (41) 

the kinetic part T 2 represents the efforts of the traders to change prices [7J. The 
intensive exchange of information in the world of finances is one of the main sources 
determining dynamics of prices. The potential part V(i?, t) of H describes the interac- 
tions between traders as well as external economic conditions and even meteorological 
conditions [7]. 

In the particular case 

H = ±f 2 

the solution satisfying the initial condition $ (i^O) — P (im) ^ s 

*(i§o.*)= e " i * f V(^o)- 

From the relation 

e-**% = t F[<p) ( T | g ) e -^ f2 $ fe 

fe=-<2 

= E F[ip] ( T g 5 ) e-^TOTK7$ fc 



obtained from ([33| and (f34|) . it follows that 



■ t ifi2 _ - t+fOOOOTTM ™2 

C 2 ^ Iy5 = C 2 f </3 

whence 

= (^,e"^ f V) =$(^o,t + 400007r M ) 

that is, the time evolution is periodic with a period of 400007T/1 seconds. This mathe- 
matical result, which does not correspond to a realistic case, shows that V(R, t) plays a 
key role in the description of the price evolution. 

The total effect of market information affecting the stock price at a certain time 
determines either the stock price's rise or the stock price's decline. In order to illustrate 
our model, by following [8] , we consider an idealized model in which we assume two type 
of information appear periodically. We choose a cosine function cos ut to simulate the 
fluctuation of the information and use the Hamiltonian 

H = ±f 2 +PRcoscjt. (42) 

where /3 is a constant. The solution of the Schrodinger equation 



<9$ 
l dt 



1 

— T + PR cos tot 
2/i 



$ (43) 
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^1 



\<S>(x,t)\ 2 



l*OM)| 2 



■ -.ill 



100 100 



100 100 



100 100 



l*0M)l 2 



, ..Tt..1 



100 100 100 100 



I4-0M)! 2 



100 100 



Figure 2: The probability |$ ( ^ , t) | 2 to have a rate of return equal to n% for t — 0, 1800, 
3600, 7200, 14400 and 28800 in the particular case a = 0.2, /m=1, /3 = 1/10, w= 1/10000. 



can be obtained by using a program in Mathematica (see the Appendix). We assume 
that at the opening time (t = 0) of the stock market, the wave function describing the 
rate of return is the function "f a (corresponding to a certain a), that is, 

$ (TM> )=7 Q (TM)- (44) 

The role played by the kinetic part of the Hamiltonian and the role played by the market 
information depend on the values of the parameters /i, /3 and uj. The distribution of the 
probabilities | SI/ (rz, t) | 2 corresponding to the possible values of the rate of return at 
certain moments of time are presented in the figure 2. We start from a stock market in 
equilibrium with the rate of return described by the Gaussian 70.2- After 1800 seconds 
the most probable rate of return is —2% and after 3600 seconds —5% and —4%. The 
rate of return comes back to equilibrium after 7200 seconds. After 14400 seconds the 
values —3% and —4% are more probable than any other value. After 28800 seconds the 
most probable rate of return is 4%. 



6 Conclusion 

Our finite-dimensional version, based on the mathematical formalism used in the case 
of finite quantum systems, keeps the essential characteristics of the model proposed by 
Zhang and Huang in [8 . It takes into consideration the discrete nature of the quantities 
used in the field of finance, and the mathematical formalism is simpler: 

- we use finite sums instead of integrals 

- in our approach any non-null function can be normalized 

- the spectra of our operators are finite sets 
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- our operators are defined on the whole Hilbert space 

- the notion of Hermitian operator can be defined in a simpler way 

- the Schrodinger type equation we use is equivalent to a system of d first order 
differential linear equations with d unknown functions. The Schrodinger equation 
used in [8] is a second order differential equation with partial derivatives. 

Appendix 

The time evolution of the probabilities t)\ 2 corresponding to the possible values 
of the rate of return can be obtained by using the following program in Mathematical 

time = 3600; alpha = 0.2; mu = 1; beta = 1/10; omega = 1/10000; d = 21; q := (d - l)/2 
g[n_] := N[ Sum[Exp[-Pi alpha (m d + n)"2/d], {m, -Infinity, Infinity}]] 
gamma := Normalize [Table [ g[n], {n, -q, q}]] 

T[n_, m_] := N[(l/(100 d) ) Sum[k Exp [2 Pi I (m - n) k/d] , {k, -q, q>] ] 
H[n_, m_, t_] := (1/(2 mu) ) Sum[ T[n, k] T[k, m] , {k, -q, q}] 

+ beta (n/100) Cos [omega t] DiscreteDelta [n - m] 
eqns = {Table[ phi[n]'[t] = Sum[- I H[n, m, t] phi [m] [t] , {m, 1, d}] , {n, 1, d}] , 

Table [phi [n] [0] == gamma [ [n] ] , {n, 1, d}] } 
ndsolve = NDSolve [eqns , Table [phi [n] , {n, d>] , {t , 30000}] 

Phi = Normalize[ Table [Evaluate [Abs [phi [n] [time] ] /. ndsolve] [ [1] ] , {n, 1, d}] ] 
ListPlot [Table [{n, Phi [ [n + q + 1]]"2}, {n, -q, q}] , Filling -> Axis] 

The reader can easily change the values of parameters time, alpha , mu, beta and the 
function Cos [omega t] in order to investigate other cases. 
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